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I. INTRODUCTION 



The B ->■ K^£ + £~ decays play a very important role in heavy flavor physics. At the 
quark level, these decays involve the flavor- changing neutral current (FCNC) of the b — >■ s 
transition, which is a purely quantum loop- mediated effect in the Standard Model (SM). 
Therefore, these decay modes have been proposed to test the SM predictions In addition 
to the branching ratio, several observables of the B — > K*£ + £~ decay, such as the longitudinal 
polarization fraction, the forward-backward asymmetry {Aps), the isospin symmetry, and 
the transverse asymmetry, have been proposed to probe possible new physics (NP) 
Various NP models thus have been scrutinized for their effects on these observables |3j. 

A few years ago, the forward-backward asymmetry of B — > K*£ + £~ was first observed 
by the Belle Collaboration J4| . The BaBar Collaboration also published its results in this 
channel earlier this year 0,JJ. Recently, the Belle Collaboration updated its measurements 
in B — > K^£ + £~ decays (7J. In these experiments, the forward-backward asymmetry is 
measured as a function of q 2 = Mfic 2 , the invariant mass of the lepton pair. In comparison, 
BaBar only has two q 2 bins of data while Belle has six. Their fitted Afb spectrum is generally 
higher than the SM expectation in all q 2 bins. This inspires us to do more investigations on 
these decays and see whether some NP model can better explain the experimental data. 

In this paper, we consider a class of family non-universal Z' models that induce FCNC's 
at tree level [8|. In such models, fermions in different families have different couplings to the 
Z' boson in the gauge basis. After rotating to the physical basis, off-diagonal couplings are 
generally produced, inducing FCNC's at tree level. These FCNC couplings are subject to 
strong constraints from low-energy experiments. Phenomenological aspects of such models 
have been extensively analyzed by various groups in recent years (9-13]. In particular, the 
possible Z'-b-s coupling has received a lot of attention because it may explain some of the 
puzzling B physics data. Based on the previous analysis, we study whether the recently 
observed B ->■ K^£ + £~ data can be accommodated within this model as well. 

This paper is organized as follows. In Sec. [Til we first review the B K ( -*H + £~ de- 
cays in the SM. In the course, we define quantities relevant for the calculations, such as 
form factors, effective Hamiltonian, explicit formulas of the amplitudes, decays widths, and 
forward-backward asymmetries. In Sec. lIIIl we describe the Z' model with tree-level FCNC's 
and deduce its effects on the B — > K^£ + £~ decays. We then use the observables to con- 
strain the model parameters. We find that the observed data in B — > K^*H + £~ can be 
accommodated in such a Z' model. We also predict the range of Br(B s — > based on 

the constrained parameter space. Finally, we summarize our findings in Sec. IIVI 
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II. B -> K(*h + £- DECAYS IN THE STANDARD MODEL 



A. Parametrization of the hadronic transitional matrix elements 



For the semileptonic decays investigated here, they involve hadronic matrix elements 
representing the B — > transitions. Therefore, we first define the B — >■ K form factors 
as follows: 



{K(p)\s^b\B(p B )) = f + {q 2 )\{pB+p) 



2 2 

m B — m K 



(K(p)\sa^q u b\B(p B )) = i{p B + p) p q 2 - q p {m 2 B - m 2 K ) 



f T (q 2 ) 



B 



m 



K 



h{q 2 )q» 



(i) 



fn B + m K 

where q = p B — p is the momentum transfer to the lepton pairs. The B — > K* transitional 
form factors are defined as: 

2V(q 2 ^ 



(K*( Pl e*)\qYb\B(p B )} 
(K*(p,e*)\q^ 5 b\B(p B )) 



(K*(p,e*)\qcr^q u b\B(p B )) 
(K*(p,e*)\qa^ l5 q u b\B(p B )} 



m B + rr%K* 



e^ pa e: PBpPa , 



2im K ,A Q {q 2 ) e —^-q^ + i(m B + m K *)Ai(q 2 ) 



-iA 2 (q 2 



e* ■ q 



q 



(Pb + P) 



M _ m B 



m 



K* u 



m B + rriK* 
-2iT x {q 2 Y vp °elp Bp p^ 

T2(q 2 ) [(m 2 B - m 2 K *)e*» - (6* ■ q)(p B +pY] 



+T 3 (q 2 )(e* ■ q) 



2 2 

m B — m K 



-{pb+pY 



(2) 



In the calculations of the semileptonic decays, we need the q 2 dependence in the form factors. 
For B — > K* transitions, we adopt the dipole model parametrization for the form factors: 

F(0) 



F(q 2 



a(q 2 /m B ) + b(q 2 /m B ) 2 



(3) 



where a and b are parameters to be determined. We calculate the form factors in the PQCD 
approach [14| near the q 2 = region, where the K* meson recoils very fast, and determine 
their values at some points. Then we extrapolate our results to the entire kinematic regime 
through fitting. Our results in the PQCD approach as well as those obtained using the QCD 
sum rules (QCDSR) 15| are listed in Table [H In our calculations, we will mainly use the 
PQCD results. The QCDSR results are included only as a comparison because we do not 
have the explicit errors on the QCDSR results . 

For the form factors of B — > K transition, we adopt a different parametrization: 



F(q 2 ) = F(0) exp [ Cl (q 2 /m 2 B ) + c 2 (g 2 /m|) 2 + c 3 (q 2 /m 2 B ) 3 } 



(4) 
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TABLE I: B ->■ K* form factors in PQCD approach and QCD sum rules (QCDSR). 





PQCD QCDSR [15] 




PQCD QCDSR [15] 


V(0) 


0.26 


0.458 


71(0) 


0.23 


0.379 


a(V) 


1.75 


1.55 


a(Ti) 


1.70 


1.59 


b(V) 


0.68 


0.575 


b(Ti) 


0.63 


0.615 


A o (0) 


0.30 


0.470 


T 2 (0) 


0.23 


0.379 


a{A ) 


1.72 


1.55 


a(T 2 ) 


0.71 


0.49 


b(A ) 


0.62 


0.680 


b(T 2 ) 


-0.19 


-0.241 


MO) 


0.19 


0.337 


Ts(0) 


0.20 


0.261 




0.79 


0.60 


a{T 3 ) 


1.58 


1.20 


KM) 


-0.09 


-0.023 


KT 3 ) 


0.49 


0.098 


A 2 (0) 




0.283 








a(A 2 ) 




1.18 








b(A 2 ) 




0.281 









TABLE II: B — > K form factors in light cone sum rules with the parametrization, Eq. @. 





F(0) ci 


c 2 c 3 


U(Q 2 ) 


0.319 1.465 


0.372 0.782 


Hi 2 ) 


0.319 0.633 


-0.095 0.591 


Mq 2 ) 


0.355 1.478 


0.373 0.700 



because the authors of Ref. find that in their fitting the extrapolation of the dipole 
parametrization to maximum q 2 is prone to reach a serious singularity below the physical 
cut starting at q 2 = m 2 B . The values of the parameters in the B — > K form factors 16( are 
listed in Table [III 



B. Effective Hamiltonian and decay amplitudes 

At the quark level, the B — > K(*H + £~ decays are dominated by the b — > si + i~ transition, 
the Hamiltonian for which is given by 

c 10 

n cS = --/§ CiMOiO*) , (5) 

* 1=1 
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where Vtb and Vt s are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and Ci{[x) 
are the Wilson coefficients evaluated at the scale /i. The local operators 0«(/i) are given by 



17| 



01 = (s a C a ) V - A (Cf3bp) V -A, 2 = (s a cp) V -A(cpb a ) v _ A) 

03 = (sgb a )v-A y](qpqp)y-A, 04 = {s a bp) V -A y^\[qpq a )v-A, 



Q 1 



05 = (s a b a ) v ^ A ^2(qpqp) v+A , Q 6 = (s a bp)y-A^2(qpqa)v+A, 



q q 
°7 = ?r^(l + 75)^ + ^^(1-75)^- 



9 = -fL(£^£)(s^(l - 7fe)6), 1O = ^(£7m75^)(^(1 - 7*)&) , (6) 

Z7T Z7T 

where a and /3 are color indices, q = u,d, s,c, {q^v-AiMilv-A = [qij^ -75)92] [937j*( 1- 
j)q 4 ), and (q^v-Afaq^v+A = [9i7 M (l - 75)^2] [g37 M ( 1 +75)94]- 

With the above Hamiltonian, the amplitude of b — > s£ + £~ transition can be written as 



+ / 



em V tb V*{cf (9 2 )[^(1 - 75)&][^] + 0io[^(l - 7fe)&][*/Vl 



2\/2 vr 

-2m 6 0f [sz<v^(l + 75)6] \hn\ - 2m s Cf [sia^{\ - Ts )6] [£7^] } , (7) 

where mj is the 6 quark mass in the MS scheme. The Wilson coefficients 0| = 07—05/3—06 
and 0g contain both the long-distance and short- distance contributions: 

9 C V) = CM + Y pert (g 2 ) + Y LD (g 2 ) . (8) 

Here Y pert represents the perturbative contribution, and Yld is the long-distance part con- 
taining contributions from the resonant states and can be excluded by experimental analysis. 
Thus we will not include Yld i n our calculation, and 

9 C V) = 9 (/i) + Y pcrt (g 2 ) , (9) 

with the detailed form of Yp ert given in Ref. 



The B — > K*£ + £ decay is more complicated because of its polarization structures in the 
final state. We will use the helicity basis. By re-expressing the metric tensor 

^ = -^ e ,(A)e:(A) + M^, (10) 

A 

we can decompose the amplitude A{B — > K*£ + £~) into Lorentz- invariant leptonic part 
L(L/R,X) and hadronic part H(L/ R, A): 

A(B^K*£+r) = L^L)H v (L)g^ + L^R)H v (R)g^ 

= - J2 X)H(L, L ^ X ) R ( R i A ) • ( n ) 
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The details have been given in Appendix C of Ref . 19[ . The explicit formulas of the functions 
L(L/R, A) and H(L/R, A) are listed in Appendix IA1 



C. The decay widths and branching ratios 



With the form factors given in Sec. Ill Al and Eq. ([7]), we obtain the dilepton spectrum of 
B -> K£ + £~ as 



dTi(B 

dq 2 



G\ | Vtb 1 2 1 Vt* \ 2 cv 2 ™ A 3 / 2 



is I u em' 



{\Ci f+(q 



2\|2 



+ 



^ ff / + ( g 2 ) + 2C7CffK + mS 



(12) 



where 



A = {m 2 K * + m| - g 



2\2 



Am 2 B m 2 K * 



[m 



B 



III 



K* 



q 2 ) 2 - 4m 2 K *q 2 . 



(13) 



For the B — > K*i + l~ decay, we define the direction opposite to the momentum of K* 
meson in the rest frame of the B meson as the +z direction. In the center-of-mass (CM) 
frame of £ + £~, 9\ is defined as the angle between the z axis and the momentum of £~ . In 
the experiment, the K* meson usually decays to the Ktt final state. We define the angle 
between the decay plane K* — > Kir and the plane determined by £ + £~ as <fi. Combining the 
leptonic amplitudes, the hadronic amplitudes, and the phase space all together, the partial 
decay width of B — > K*£ + £~ is given by 



dTi(B -»■ K*> 

vg 

10247r 4 m| 



10247r 4 m| 



d cos exd(j)dq 2 \Ai{B K*( 



d cos M(M? 2 (|L(L, i)H(L, i)\ 2 + \L(R, i)H(R, i) 



(14) 



where i = 0, + or — denotes the three different polarizations of the K*. 

After integrating out 9\ and <fi in Eq. ( Tl4l) . one obtains the dilepton spectrum of B 
K*£ + £~ decay as: 



dTi(B -»■ 



dq 2 



V\q 2 
967r 3 m 3 3 



|if(L,z)| 2 + |F(it;,z) 



(15) 



In Sec. Ill B\ one can find that among the Wilson coefficients only Cg has the q 2 depen- 
dence. The dilepton spectra of B — > K^£ + £~ decays are shown in Fig. HJ with and without 
^pert(? 2 ) in Cg S being included. After further integrating out the q 2 dependence, we obtain 
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FIG. 1: ^-dependence of the branching ratios of B — > K£ + £~~ (left plot) and B — > K*£ + £~ (right 
plot) decays. In the left plot, the red solid (blue dashed) curve stands for the dilepton spectrum 
with (without) the Y pevt (q 2 ) part included in Cf. The right plot shows the spectrum predicted 
in PQCD and QCDSR with and without the Y pert (q 2 ) part in Eq. Q. The black solid (red long 
dashed) curve is the PQCD results with (without) Y pert (q 2 ) and the blue short dashed (pink dotted) 
curve is the QCDSR result with (without) Y per t(q 2 )- In the curves where Y per t(q 2 ) is included, a 
kink shows up because it is a piecewise function. 



the total branching ratios: 



Br(B KIH-) = i (4 - 7 °-" 7?) X 1(r? Part " ° f inClUd6d) ' (16) 
1 (4.45tJ;g) x 10" 7 (q 2 part in Cf excluded), 

Br(B = { (16 - 5 "" ) X 10 ~ 7 ^ in ° f included )' (17) 



'15.8±H) x KT 7 {q 2 part in Cf excluded). 
These predictions are to be compared with the experimental results Q: 

Br(B -> K£+r) = (4.8^;1 ± 0.3) x KT 7 , 

Br(B -»■ K*£ + r) = (10.7ll;5 ± 0.9) x 10~ 7 . (18) 

From Fig.lU Eq. (fTBl and Eq. ( JTTjl . one finds that the Y peTt (q 2 ) piece in Cf has a small effect 
on the branching ratios in comparison with other uncertainties. To simplify the notation, 
we define C g = Y pcrt (q 2 ), and thus Cf = Cg + C' 9 . The differential branching ratio of 
B — > Ki + i~ is then decomposed into the following form 

dBr(B = ^ B ' 1 + |c<9Cff|2jB2 + |c?ff + ^[cf cf*]^ 

= |C 10 | 2 ^ + [|C 9 | 2 + |C^| 2 + 2Re[C 9 C' 9 *]]B' 2 + \Cf\ 2 B! 3 + 2Re[(C 9 + C£)Cf .(19) 
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After the integration over q 2 , Eq. (fT9l can be rearranged as 

Br(B K£ + £~) = \C W \ 2 B 1 + \C 9 \ 2 B 2 + \Cf \ 2 B 3 + 2Re[C 9 Cf*]B 4 + 2Re[C 9 ]B 5 

+2Re[Cf]B 6 + B 7 , (20) 

where B 5 (B 6 , B 7 ) contains the integration of Re[C' 9 ]B' 2 {Re[C' 9 ]B' A , \Cg\ 2 B' 2 ). Similarly, 
Br{B — > K*£ + £~) is decomposed as 

Br[B K*£ + t) = \C 10 \ 2 Bl + \C 9 \ 2 B* 2 + \Cf\ 2 B* + 2Re[C 9 Cf*]Bl + 2Re[C 9 ]B* 

+2Re[cf ]b; + b; . (21) 

The values of B [ * ] with j = 1, 2, 3, 7 are, in units of 10" 8 (1(T 7 ), 



B 1 


_ i 9O+0.30 
— l.ZO_ 23 ■ 


B 2 


= B X 


, B 3 = 


4 41 + 1 - 44 


B,= 


9 qq+0.71 
z -°°-0.39 > 


B 5 


_ n oi +0.09 
— u.oi_ 005 , 


B 6 


= 0.5£ 


>+0.19 
'-0.10 ) 


B 7 = 0.18; 


fO.04 
-0.03 ) 




B{ 


- n4i+ - 20 


B* 2 


= B* 


, B* 3 = 


1 9 74+ 6 - 35 

J-^- '^-4.86 


TD* 




B* h 


- o nq+ - 04 


b; 


= O.U 


j+0.10 
5 -0.10 ) 


B* = 0.04 


+0.02 
-0.02 • 





(22) 

These values will be used to constrain the couplings in the Z' model later. From Fig. [1] and 
Eqs. (TAT]) to flAT2|) . one can find a pole at q 2 = in d J Br( J B K*£ + £~)/dq 2 . That is why 
is much larger than the others. 

D. The forward-backward asymmetry 

The differential forward-backward asymmetry of B — > K*£ + £~ is defined by 
dAps f 1 d 2 T f° d 2 T 



dq 2 J dq 2 d cos 9i J_ 1 dq 2 d cos 6^ 

while the normalized differential forward-backward asymmetry is defined by 



dA FB 


dA FB 
dq 2 


dq 2 


dF 




dq 2 


3 





-\H(L, +)\ 2 + \H(R, +)| 2 + \H(L, -)\ 2 - \H(R, -)\ 2 

4 \H(L, 0)| 2 + \H(R, 0)| 2 + \H(L, +)| 2 + \H(R, +)| 2 + \H(L, -)| 2 + |JT(i2, -)| 2 



(24) 



Substituting the expressions in Eqs. (1A7|) to (1A12I) into Eq. (1241) . we get the explicit expres 

dA FB 3N(q 2 ) 



sion for d d ^ 2 B as follows: 



dq 2 4D(q 



2^ > (25) 
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where 



N(q 2 ) = \V tb \ 2 \V;fG 2 F a 2 m v^q 2 {-Re[C 10 )Cfm b [(m B 

+ {m B - m K *)T 2 (q 2 )V(q 2 )] + Re[CfC* w ] [-q 2 V(q 2 )A 1 (q 2 )] } , 

D(q 2 ) = 2n 2 {q 2 f [\H{LM 2 + \H{RM 2 + \H{L,+)\ 2 + \H{R,+)\ 2 

+ \H(L,-)\ 2 + \H(R,-)\ 2 ] . (26) 

In the above expression, terms suppressed by m s are dropped for simplicity. As can be 
explicitly checked, the pole in the dilepton spectrum at q 2 = disappears in the denominator. 




FIG. 2: The forward-backward asymmetry for B — > K*£ + £ with form factors given by the PQCD 
approach. The black solid curve is given with the SM CV and the red dashed curve is given with 
C 7 = 0. 



According to Eq. (126 p . the numerator of dAps/dq 2 is zero at q 2 = because of the 
common factor q 2 , while the denominator has a non-zero value because its common factor 
(q 2 ) 2 cancels with the (q 2 ) 2 factor arising from Eqs. ( lA8l) . ( TA~9l . (lAllj) and (1A12I) . Thus 
dA FB /dq 2 = at q 2 = 0. In the SM, Cf < 0, Cf > 0, and C w < 0; thus the first term 
in the curly bracket of N(q 2 ) is negative and the second term is positive. In the regime 
where q 2 is near zero, the first term gives the dominant contribution since the second term is 
suppressed by the small q 2 . Therefore, the sign of dAps/dq 2 is determined by the first term 
and gives a negative value. As q 2 increases, the second term becomes dominant. There exists 
a point where dAps/dq 2 becomes zero, the so-called forward-backward asymmetry zero. The 
position of the zero is determined by C% s and Cg ff , for the form- factor dependence drops at 
the leading order As q 2 becomes even larger, the effect of the overall factor y/\ becomes 
crucial. Eq. ( fl3|) tells us that A = at the largest recoil where q 2 = (mB — mK*) 2 - Therefore, 
dAps/dq 2 falls back to zero at the end of the kinematic regime. All these behaviors of 
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dApB/dq 2 can be observed in Fig. [2j The red dashed curve is drawn with the contribution 
of only the second term in the curly bracket of N(q 2 ). It shows the importance of C% in 
the low q 2 regime. 

However, the latest Belle data 0] do not show an obvious zero for dAps/dq 2 , and the 
values at all q 2 are consistently higher than the SM expectation. A common solution is to 
flip the sign of Cj S as it is still consistent with the constraint from B — > X s ^ data. In the 
next section, we offer an alternative solution in the family non-universal Z' model. 



III. CONSTRAINTS ON THE COUPLINGS IN Z' PHYSICS 



A. b -»• s£ + £~ in the Z' FCNC model 



In the appropriate gauge basis, the U(l)' currents are 

3% =g'J2 ^\4 l Pl + 4 r Pr)^ (27) 

i 

where i is the family index and ip labels the fermions (up- or down-type quarks, or charged 
or neutral leptons), and Pl,r = (1 -F 7s)/2- According to some string construction or GUT 
models such as Eq, it is possible to have family non- universal Z' couplings. That is, even 
though ef' R are diagonal, the couplings are not family universal. After rotating to the 
physical basis, FCNCs generally appear at tree level in both LH and RH sectors. Explicitly, 

B** = V^Vl, B** = V^Vl. (28) 

Moreover, these couplings may contain CP-violating phases beyond that of the SM. 
In particular, Z'bs couplings can be generated: 

£|cnc = -g'{B^ h s Llll b L + B? b s Rllx b R )Z'» + h.c. . (29) 

The couplings in Eq. ( 129]) lead to extra contributions to the b — > s£ + £~ decay at tree level, 
mediated by a virtual Z' boson. The amplitude is given by 

g' 2 



M% 



[B L sh s Llli b L + Bf b s Rl ,b R ) [B L a i Ll n L + B*I Rl n R ) . (30) 

There are thus four types of operators, Oll, Olr, R l, and RR . The above amplitude can 
be derived from an effective Hamiltonian 

w£ = ^(ffe^br, + pf b s Rl ,bMhrh + pgtxfiR) , (si) 



where 



L,R _ 9 Mz D L,R ( on\ 

Pff ' ~ gM z ff ' [ ' 
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TABLE III: Values of Wilson coefficients Ci(rrib) in the leading logarithmic approximation, with 
mw = 80.4GeV, [i = rrn, jP0 i e [17]. 



Ci 


c 2 


c 3 


c 4 


c 5 


c 6 




C10 


1. 107 


-0.248 


-0.011 


-0.026 


-0.007 


-0.031 


-0.313 4.344 


-4.669 



and g is the coupling associated with the SU(2)l group in the SM. Throughout this analysis, 
we ignore the renormalization group running effects due to these new contributions because 
they are expected to be small. 

B. Constraints from the B — > K^*H + i~ decays 

For the purpose of illustration and to avoid too many free parameters, we assume that the 
FCNC couplings of the Z' and quarks only occur in the left-handed (LH) sector. Therefore, 
pf b = 0, and the effects of the Z' FCNC currents simply modify the Wilson coefficients Cg 
and C10 in Eq. fl5]). We denote these two modified Wilson coefficients by Cg S ' z and C^,', 
respectively. More explicitly, 

Re[Cf*'] = Re[Cf]-^ Re[p " Ml + P ^ 



V tb V t * s a em 

Im[Cf> z ) = Irn[Cf]-^ Im[p " M + p ^ 



Re[Cf ] = C 
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V tb V t * s a em 
V tb V t * s a em 



v tb v ts ot em 

For simplicity, we further assume that p^ b is real. Then the imaginary part of Cg S will not 
be affected by the Z' model, and Cf§ is still a real number. 

First, we consider the constraint from the spectrum of dApB/dq 2 . In order to fit the 
experimental data, a sign flip is needed for dAps/dq 2 near the q 2 = regime. People 
usually consider the flipped-sign solution with CV = — Cf M , because it is still allowed by 
the B — > Xgj data. However, an alternative solution is to flip the signs of C| and C\o 
instead, as is possible in our model. Below Eq. fl26|) . it is noted that in this regime the term 
proportional to Re[Cio]C^ E dominates. Therefore, one can flip the sign of Ciq: 

Re[C?o] > . (34) 
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Moreover, in order to keep the second term in the curly bracket of N(q 2 ) to have the correct 
behavior, we also need to flip the sign of Re[Cg S ]. Thus, we require 

Re[Cf' z ) < . (35) 

Eqs. (13~4"|) and (I3"5"j) are the constraints from the dAps/dq 2 spectrum obtained by the Belle 
Collaboration (see Fig. 1 in Ref. 0]). 

Next, we consider the constraints from the branching ratios of B — > K^£ + £~ decays. 
These constraints are obtained in the following way. After including the contributions of 
Z', the upper (lower) bound of the theoretical predictions should be greater (smaller) than 
the experimental lower (upper) bound at the 2a level. When we deal with the experimental 
data, we add the statistical and systematic errors in quadrature. With Eqs. f f2T]) and fl22l) . 
we have the following branching-ratio constraints: 

aff(|C£| 2 + \C* | 2 ) + B^\Cf\ 2 + B£Re[C?Cf*\ + 2&>Jfe[C7f ] 

+B^Re[Cf } + B£ > Brg p - 2a« , (36) 

Bu(\ c w\ 2 + \C* | 2 ) + B$\Cf\ 2 + B$Re[C?Cf*\ + B^Re[C 9 z ' ] 

+B$Re[Cf] + flW < Sr£J + 2a« , (37) 

where quantities with a star in the superscript are for the B — > K*£ + £~ decay, and the letters 
"it" and "Z" in the subscript represent the 1-a upper and lower bounds of the correspond- 
ing quantity B^*\ respectively, and Br&q, denote the central values of the B — > K^£ + £~ 
branching ratios. 

Moreover, Cg = Cg + x with Y peTt (q 2 ) excluded, and C^' = C w + y, where 

vtbv ts a em 

v = ^Mp L sM - Pu) m) 

U V tb V t * s a em ' 1 } 

Then Eqs. (I3"6"j) and fl3T|) can be rearranged as 

(x + T«) 2 + B« (y + C 10 ) 2 + C« > BrW - 2<r, w , (40) 
+ T/* } ) 2 + B«(y + Cio) 2 + C<* > < Br« + 2a« , (41) 



where 



W . . 2 <A + B$ fl Cf + flW 



(*) _ r(*) r<2 | d(*) //-yeff\2 , r?(*) r< , r>(*) , r>0) 



r<(*) _ rW n 2 _L_ Rv ' //mett\2 , p>\*) syctt /~t , R l*J <o , R (.*J <o 
~~ 9 + jW 7«A ( -°7 J + i V/Z°7 5u/l & 6ii/Z 7 

)(*) _ r(*) /T-iC*) 

J 7«/i lu/ZV u/Z 



( T W) 2 . (42) 
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FIG. 3: The constraints from branching ratios of B — > K^i + £~ decays. The areas outside the 
red solid and yellow short dashed circles are determined by Eqs. (|45p and (|47p . respectively. The 
areas inside the pink and blue circles are determined by Eqs. (|46|) and (|48|) . respectively. The areas 
to the left of the line x = —Cg and above the line y = — C\q are determined by Eqs. (|43[) and 
(|44p . respectively. The black dot is where both Cg and Cio flip signs from their SM values. The 
two rectangles, corresponding to SI (the large rectangle) and S2 (the small rectangle) in Case III, 



are the constraints given by Ref. [I3j|. One can see that their constraints are consistent with our 
constraints from the branching ratios. However, their constraints are not enough to change the 
signs of Cg and C\q. 



Substituting all the numerical values in Eqs. (l3lj) . fl35|) . ( 140)) . and fl4T]) . we have 

x < -4.344 , (43) 

y > 4.669 , (44) 

1.58(x + 3.99) 2 + 1.58(y - 4.669) 2 - 37.88 > , (45) 

1.05(x + 4.01) 2 + 1.05(y - 4.669) 2 - 59.58 < , (46) 

0.61{x + 3.89) 2 + 0M(y - 4.669) 2 - 6.38 > , (47) 

0.26(x + 4.11) 2 + 0.26(?/-4.669) 2 - 12.81 < 0. (48) 

Eqs. ( H3"j) -( H8"|) give the constraints on x and y, which are shown in Fig. [3J The common area 
of the above six conditions is outside the red solid circle and inside the blue long dashed 
circle, to the left of the solid vertical line x = —Cg and above the solid horizontal line 
y = —Ciq. This area gives 
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(Br*+2a* u -Cf)/Bl 



-C w < y < ^(Br* xp + 2a* u - Cfi/B^ - C 10 
With Eqs. (155]) and (135]), we have 



(49) 



(Br* + 2a* u ~ Cf)/B* u - C w - C 9 K< Retfdpg 



< 



(Br* exp + 2< - C?)/B* u - T? - C 
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K 



[C 10 - C 9 ] K < Re[ P L s Mi 



< 



cxp 



2a* u -Ct)/B* l -T l * + C 



in 



(50) 



with /C = (VJb^*a em )/(47r). In the quark sector, the coup ling s in Eq. fl29l also lead to a 
NP contribution to B® — B® mixing at tree level. In Refs. [20|, |2l|, it is assumed that only 
the LH sector of quarks has family non-universal U(l)' couplings, as in the current analysis. 
Thus, only the LH interaction in Eq. ( |29|) contributes to the B® — B° s mixing. They find 
that one can reproduce the measured value of AM S if 

Psb Z 1(T 3 . (51) 
As a rough estimate, here we take p^ b = 1CT 3 . Together with Eqs. ( )38|) . (|39|) . and 



Vtb = 0.999176 , V u 



-0.03972 , a em 



1/137 



we obtain 



0.27 < p\ { < -0.11 
-0.08 < p% < 0.09 . 



(52) 



(53) 
(54) 



We should emphasize that these parameter ranges are obtained with some assumptions and 
the current data. In particular, we have used a particular value of p L sh for our illustration. 
Once new experimental data or theoretical inputs are available, these constraints can be 
easily updated with our formulas. In Figj3]we also give the constraints from Ref. 13]. In 
their paper the authors gain the constraints by making the experimental and theoretical 
values of B — > X s l + l~ agree with each other in la. However, we get our constraints in 2a. 
For a comparison, in Fig|3]we simply extrapolate their results to 2a. One can find that if 
we drop the constraint conditions of flipping the signs of Cg and Cio, we agree with each 
other. However, if the Afr is expected to behave as how we get constraints ( 1341) and ( l35l) . 
the constraints in Ref. 



13| are too tight to satisfy the conditions. 
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0.8 
0.6 
0.4 
0.2 


-0.2 

5 10 15 

q 2 

FIG. 4: Forward-backward asymmetry in QCDSR (red dotted line) and PQCD (black solid line) 
with Cg and Cio flipping their signs (the black dot in Fig. ED. The points with error bars are the 
experimental results from the Belle Collaboration Qj. 

In Fig. HI we use the black dot from Fig. |3l where both C| ff (g 2 ) and C\q flip signs from 
their SM values, to predict the dAps/dq 2 spectrum in our model. Since Cg S (q 2 ) is q 2 - 
dependent, the plot in Fig. [3]is plotted with Cg and C w flipping their signs. The points that 
flip the signs of Cg S (q 2 ) and C\q should be very close to this point. It is interesting to note 
that the red dotted curve in Fig. [4] is identical to the usual flipped-sign solution. This is not 
surprising because flipping the signs of both Cg S (q 2 ) and C w simultaneously is equivalent to 
flipping the sign of Cj S , which can be seen from Eq. (1261) . This indicates that by considering 
only the branching ratios and forward-backward asymmetry of the B — > K^*H + i~ decays, 
it is insufficient to determine which operators are significantly modified by the NP. 

Now a comment on the form factors is in order. Because of the nonperturbative effects, we 
cannot get good results for the form factors when q 2 is large. In either PQCD or light cone 
sum rules, the form factors are obtained in a region where q 2 is small and then extrapolated 
to the entire kinematical region through fitting. As a result, it is a question whether the 
form factors can be described well by the parametrization formula in the large q 2 region. In 
fact, the accuracy of the parametrization formula becomes worse as q 2 increases. Therefore, 
we do not think the theoretical predictions at large q 2 are reliable enough. This may explain 
why the experimental values are still a little larger than the theoretical predictions in the 
large q 2 regime, as shown in Fig. |H 

A closely related decay mode to the current analysis is the B s — > fi + fi~ decay. This mode 
has been searched for with great interest at Tevatron. The upper bounds on the branching 
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ratio at 95% confidence level are given by its two experimental groups as 



Br(B a -> < 5.8 x lCT 8 (CDF) , 



Br(B s < 1.2 x 1(T 7 (D0) [23] 



(55) 



The branching ratio of B s — > fi + fi is affected in our model. With the inclusion of the Z' 
contribution, the branching ratio is given by 171 ] 



Br(B s — > 



T Bs^:fB s <m Bs Jl 



4m?, 



47T 



71} 



B 3 



F f 2 |T/*T/ , 

V tb V ts | 

L 



a 



2n sin 9 W 



w 



As(p l 



V&Vu 



(56) 



where all the functions and symbols are defined in Ref. 17] . With the constraints in Eq. (149|) . 
we find that the upper bound for this branching ratio is 



Br(B a -> ^i + /i^) < 7.9 x 10^ 9 . 



(57) 



Note that the upper bound of the range is still smaller than the current upper bound given 
by CDF Collaboration. 



IV. SUMMARY 

We have considered the contributions of family non-universal Z' models with flavor- 
changing neutral currents [Z' FCNC) at tree level in B — > K^*H + £~ decays. By requiring 
that the theoretically predicted branching ratios agree with the current experimental data 
within two a's, we obtain the constraints on the couplings in the Z' FCNC model. We 
find that within the allowed parameter space, our model has the potential to explain the 
forward-backward asymmetry of the B — > K*£ + £~ decay, as better determined by the Belle 
Collaboration recently. Moreover, our Z' model contributions flip the signs of Cg G and Cio, 
which differs from the usual new physics contributions that flip the sign of C% . Using the 
constraints, we also compute the branching ratio of the B s — > fi + n~ decay. The upper bound 
of our prediction is near the upper bound given by CDF Collaboration. 
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Appendix A: functions for the leptonic and hadronic part 



L(L,0) 


= 2^/q 2 sin^i, 


L(L,+) 


= -2v / 2\ 


' q 2 sin 2 


L(L,-) 


= -2y/2 % 


/ q 2 cos 2 


L(R, 0) 


= 


sin 9i, 


L(R,+) 


= -2v^\ 


I q A cos 


L(R, -) 


= -2y/2\ 


/ q 2 sin 2 



2 ' 
2 



its 



2 ' 
0i 
2 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 

(A6) 



H(L,0) 



iG F V tb V t * s a 

em 

8y/2~7rmK* \fq* 

+(cf - c 10 ) 



2(C 7i - C 7R )m b 



AT 3 (g 2 ) 



(m B + m K »)(m K « - m B + q )A 1 (q ) + 



(3m^+m|-g 2 )T 2 (g 2 ) 
AA 2 (g 2 ) 



(m s + m K *) 



(A7) 



em 



4\/27rg 2 
+(Cf - C 10 )g 2 



2(C 7L + C 7i? )m fe v / AT 1 (g 2 ) - 2(C 7L - C m )m b (m 2 B - m 2 K ,)T 2 (q 2 ) 



r 



[m B + m K *. 



- (m B + m K *)A 1 (q 2 ) 



(A8) 



H(L, 



H(R,0) 



iG F V tb V t * s a 

cm 

4 v / 2vrg 2 



+ (Cf - C 10 )q 2 



iG F V tb V*a 

em 

%\p2i\ViiK* ^/q 2 
+(Cf + C 10 ) 



2(C 7L + C 7R )m h V\T x {q 2 ) - 2(C 7L - C 7R )m b (m 2 B - m 2 K AT 2 (q 2 ) 
^V(q 2 ) 



(m B + m K *) 



2{C 7L - C 7R )m b 



(m B + m K *)A x {q ) 



(A9) 



XT 3 (q 2 



m B - m« K » 



(m B + m K *)(m K * - m B + q )A x (q ) 



(3m 2 K * +m 2 B - q 2 ) T 2 (q 
XA 2 (q 2 ) 



(m B + m K *) 



(A10) 



H{R, 



iG F V tb V t *a 



is "em 



4y / 27rg 2 
+ (Cf + C w )q 2 



2{C 7L + C 7/? )m 6 v / AT 1 (g 2 ) - 2(C 7L - C 7R )m b (m 2 B - m 2 K ,)T 2 (q 2 ) 



V\V(q 2 ) 
(m B + m K * 



(m B + rriK^A^q 2 ) 



(All) 
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H{R, 



iG F V tb V*a 



ts^em 



4\/27rg 2 



+(Cf + C 1Q )q 2 



2(C 7L + Cr^mbVXTtiq 2 ) - 2(C 7L - C 7R )m b {m 2 B - m 2 K ,)T 2 {q 2 



V\V(q 2 ) 
(m B + m K *) 



(m B + m K *)A 1 (q' 2 



(A12) 



[1] Jaus W and Wyler D, Phys. Rev. D, 1990, 41, 3405; Colangelo P, De Fazio F, San- 
torelli P and Scrimieri E, Phys. Rev. D, 1996, 53, 3672 [Erratum-ibid. D, 1998, 57, 3186 



arXiv:hep-ph/9510403] ; Aliev T M, Ozpineci A and Savci M, Phys. Rev. D, 1997, 56, 4260 



arXiv:hep-ph/9612480| ; Melikhov D, Nikitin N and Simula S, Phys. Lett. B, 1998, 430, 332 



|arXiv:hep-ph/9803343] . 

[2] Burdman G, Phys. Rev. D, 1995, 52, 6400 [arXiv:hep-ph/9505352| ; Burdman G, Phys. Rev. 
D, 1998, 57, 4254 [arXiv:hep-ph/9710550| ; Beneke M, Feldmann T and Seidel D, Nucl. 
Phys. B, 2001, 612, 25 |arXiv:hep-ph/0106067| ; Feldmann T and Matias J, JHEP 2003, 
0301, 074 |arXiv:hep-ph/0212158| ; Kruger F and Matias J, Phys. Rev. D, 2005, 71, 094009 
[arXiv:hep-ph/0502060] . 

[3] Hewett J L and Wells J D, Phys. Rev. D, 1997, 55, 5549 [arXiv:hep-ph/9610323l ; Ali A, 
Ball P, Handoko L T and Hiller G, Phys. Rev. D, 2000, 61, 074024 |arXiv:hep-ph/9910221| ; 
Colangelo P, De Fazio F., Ferrandes R. and Pham T. N., Phys. Rev. D, 2006, 73, 115006 
arXiv:hep-ph/0604029] ; Hovhannisyan A, HOU W S and Mahajan N, Phys. Rev. D, 2008, 



77, 014016 [arXiv:hep-ph/0701046| . 
[4] Ishikawa A et al, Phys. Rev. Lett. 2006, 96, 251801 |arXiv:hep-ex/0603018| . 
[5] Aubert B et al. [BABAR Collaboration], Phys. Rev. Lett. 2009, 102, 091803 |arXiv:0807.4TT9l 

[hep-ex]]. 

[6] Aubert B et al. [BABAR Collaboration], Phys. Rev. D, 2009, 79, 031102 |arXiv:0804.44l2l 
[hep-ex]]. 

[7] J. T. Wei et al. [BELLE Collaboration], Phys. Rev. Lett. ,2009, 103, 171801 |arXiv:0904.0770l 
[hep-ex]]. 

[8] Langacker P and Plumacher M, Phys. Rev. D, 2000, 62, 013006 [arXiv:hep-ph/0001204] . 

[9] Barger V, Chiang C W, Langacker P and Lee H S, Phys. Lett. B, 2004, 580, 186 
arX iv:hep-ph/0310073| ; Barger V, C. W. Chiang, J. Jiang and P. Langacker, Phys. Lett. B, 
2004 596, 229 [arXiv:hep-ph/0405108| ; Barger V, C. W. Chiang, P. Langacker and H. S. Lee, 
Phys. Lett. B, 2004, 598, 218 [arXiv:hep-ph/0406126] ; Arhrib A, Cheung K, C. W. Chiang 
and T. C. Yuan, Phys. Rev. D, 2006, 73, 075015 [arXiv:hep-ph/0602175j ; Cheung K, Chiang 
C W, Deshpande N G and Jiang J, Phys. Lett. B, 2007, 652, 285 larXiv:hep-ph/0604223l ; 



18 



Chiang C W, Deshpande N G and Jiang J, JHEP, 2006, 0608, 075 [arXiv:hep-ph/0606122] . 
[10] HE X G and Valencia G, Phys. Rev. D, 2004, 70, 053003 [arXiv:hep-ph/0404229] | ; Phys. Lett. 
B, 2007, 651, 135 |arXiv:hep-ph/0703270] . 



[11] Barger V, Everett L, Jiang J, Langacker P, Liu T and Wagner C, Phys. Rev. D, 2009, 80, 

055008 [arXiv:0902.4507l [hep-ph]]: la7Xiv:0906.3745l [hep-ph]. 
[12] CHEN S L and Okada N, Phys. Lett. B, 2008, 669, 34 |arXiv:0808.033T1 [hep-ph]]. 
[13] CHANG Q, LI X Q and YANG Y D, JHEP 2009, 0905, 056 [arXiv: 0903 .02751 [hep-ph]]; 

larXiv:0907.4408l [hep-ph]. 
[14] Keum Y Y, Li H n and Sanda A I, Phys. Lett. B, 2001, 504, 6 [arXiv:hep-ph/0004004] ; LU 



CD, Ukai K, YANG M Z, Phys. Rev. D, 2001, 63, 074009 e-Print: hep-ph/0004213; Keum Y 



and Li H n, Phys. Rev. D, 2001, 63, 074006 |arXiv:hep-ph/0006001] ; LU C D, YANG M Z, 



Eur. Phys. J. C, 2002, 23, 275-287 e-Print: |hep-ph/0011238 
[15] Ball P and Braun V M, Phys. Rev. D, 1998, 58, 094016 [arXiv:hep-ph/9805422| . 
[16] Ali A, Ball P, Handoko L T, and Hiller G, Phys. Rev. D, 2000, 61, 074024. 
[17] Buchalla G, Buras A J and Lautenbacher M E, Rev. Mod. Phys. 1996, 68, 1125 



|arXiv:hep-ph/9512380] . 
[18] Buras A J and Munz M, Phys. Rev. D, 1995, 52, 186 [arXiv:hep-ph/9501281| . 
[19] LI R H , LU C D and WANG W, Phys. Rev. D 2009, 79, 094024 |arXiv:0902.329T1 [hep-ph]]. 
[20] Cheung K, Chiang C W, Deshpande N G and Jiang J, Phys. Lett. B, 2007, 652,285 
[21] V. Barger, L. Everett, J. Jiang, P. Langacker, T. Liu and C. Wagner, Phys. Rev. D, 2009, 80, 

055008 [arXiv:0902.4507l [hep-ph]] : larXiv:0906.3745ft iep-ph] ; V. Barger, L. L. Everett, J. Jiang, 

P. Langacker, T. Liu and C. E. M. Wagner, JHEP, 2009, 0912, 048 |arXiv:0906.3745l [hep-ph]] 
[22] Aaltonen T et al. [CDF Collaboration], Phys. Rev. Lett. 2008, 100, 101802 [ arXiv:0712.1708l 

[hep-ex]]. 

[23] Abazov V Met al. [DO Collaboration], Phys. Rev. D 2007, 76, 092001 |arXiv:0707.3997l [hep- 
exl]. 



19 



